Mimicking Nanoribbon Behavior Using a Graphene Layer on SiC 
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We propose a natural way to create quantum-confined regions in graphene in a system that 
allows large-scale device integration. We show, using first-principles calculations, that a single 
graphene layer on a trenched region of [0001] SiC mimics i)the energy bands around the Fermi level 
and ii) the magnetic properties of free-standing graphene nanoribbons. Depending on the trench 
direction, either zigzag or armchair nanoribbons are mimicked. Fhis behavior occurs because a 
single graphene layer over a SiC surface loses the graphene-like properties, which are restored solely 
over the trenches, providing in this way a confined strip region. 

PACS numbers: 73.22.Pr,71.15.Nc,73.20.At 



Graphene has attracted enormous attention due to its 
potential for application in devices at the nanometer 
scale. The very high mobility of graphene charge carri- 
ers, exceeding 200.000 cm 2 /V.s, presents a possibility of 
fabrication of field effect transistors (FETs) as an alter- 
native to .Si-based devices [lj . As is well known, graphene 
has a linear dispersion relation with valence and conduc- 
tion band degenerate at the so called Dirac point 
This electronic behavior implies that, although pristine 
graphene can be used for applications in radio frequency 
(RF) devices [|, it has no energy gap. Consequently, 
it presents a very small current on/off ratio for transis- 
tors Thus, one of the most fundamental challenges 
for researchers in the field is to fabricate graphene-like 
structures with a band gap that i) preserve the above 
mentioned qualities; ii) the band gap is fabricated in a 
controllable way and iii) in such a way that it is possible 
to integrate them in large-scale circuits. 

A possible solution for these challenges would be to 
use graphene nanoribbons (GNR) Q. Huge efforts 
have been made with the aim of understanding the prop- 
erties of GNRs, and possible applications on electronics, 
spintronics, memory devices and sensors have been pro- 
posed in the literature j7-ll|. The gap can be controlled 
by the width of the ribbon and the high mobility is still 
preserved in these systems. Considering large-scale cir- 
cuits and the perspective of use of graphene on wafers for 
industrial applications, the process of obtaining a single 
layer of graphene via heating a hexagonal SiC surface 
until the sublimation of Si surface atoms (with tempera- 
tures over 1000 °C) l_2j is attracting attention as a viable 
option. Recently, there have been proposals for FETs 
which are based on graphene layers grown on SiC sub- 
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In this letter we propose a system that naturally cre- 
ates graphene nanoribbons on top of SiC substrates. We 
show, using ab-initio density functional theory (DFT)[lR 
17j . that a single sheet of graphene over the [0001] 
(carbon rich) surface of SiC can mimic a free-standing 
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FIG. 1. (Color online) (a) Schematic picture of a single 
graphene layer over a SiC surface containing a trench, (b-e) 
Side and top views of the local magnetization (pt(r) — P4.(r)) 
for a graphene sheet over SiC with a 20. 6 A wide trench cut 
along the T — K direction, (b-c) Anti- ferromagnetic, and (d-e) 
ferromagnetic configuration. Dark (blue) indicates excess of f 
electrons, while the light (green) indicate excess of 4- electrons. 



graphene nanoribbon (FS-GNR) in places where sub- 
strate atoms are removed, creating in this way trenches 
in the SiC 21-23}. In Fig. [TJa) a schematic picture of 
our proposal is depicted. Depending on the direction of 
the trench, either FS-GNRs with zig-zag (FS-zGNRs) or 
armchair edge shapes (FS-aGNRs) are mimicked. When 
the trench is in the T — K graphene direction, ferromag- 



FIG. 2. (Color Online) Energy bands, (a) Single sheet of graphene over a crystalline surface of SiC without a trench. Anti- 
ferromagnetic configuration: (b) FS-zGNR and (c) single sheet of graphene over SiC with a trench in the T — K direction. 
Ferromagnetic configuration: (d) FS-zGNR and (e) graphene over trenched SiC. (f-j) Self-interaction correction for cases 
reported in (a-e). The trenches considered here have a width of 20. 6 A. The solid (blue) lines for the t states, and the dashed 
(red) lines for the ! states. 



netic (metallic) and anti-ferromagnetic (semi-conductor) 
solutions are possible, exactly as one would find in the 
FS-zGNRs. For the trench along the V — M graphene 
direction, tree families of energy bands can be identified 
and associated to the 3p, 3p + 1 and 3p + 2 families of FS- 
aGNRs, where p is a positive integer 24 1. For graphene 
over the [0001] (silicon rich) SiC surface, the half-filled 
dangling-bond states pin the Fermi energy at ps 0.5eV 
above the Dirac point [25| (near the bottom of the con- 
duction band) , preventing the appearance of states that 
mimic FS-GNR near the Fermi energy, even in the pres- 
ence of a trench. 

In our geometry, a single layer of graphene is posi- 
tioned over SiC surfaces with the y/3 x -\/3-R30 ° recon- 
struction. The substrate we consider has 3 bilayers of 
SiC in the ABC stacking pattern, representing a surface- 
terminated thin slab of 4H-SiC with periodic boundary 
conditions. In this situation, there are covalent bonds 
between the graphene and the top-surface atoms that 
destroy the free-standing graphene properties (such as 
the conic dispersion around the Fermi-level), in agree- 
ment with previous ab-initio simulations [2H \2l\ as well 



as experiments performed in both C-rich [28[ and Si-rich 
surfaces j29[. We obtain the trench by removing atoms 
of the substrate from underneath the graphene layer in 
the central region, considering two possible directions, 
namely along the V — K and T — M graphene Brillouin 
zone directions. Along the V — K direction, trenches 7.5, 
11.9, 20.6 and 34.8 A wide were considered, whereas for 
the r — M direction, we consider trenches 5.2, 7.7, 10.4, 
13.0, 15.7, 18.2, 20.9, 23.5 and 26.2 A wide. 

In Fig. [5] we show the energy bands along the T — K 
direction. (The Fermi energy is shifted to 0.0 eV in all 
graphs). Fig. Ufa) represents the dispersion relation for 
the system without a trench. The two flat bands above 
and below the Fermi-level are associated with dangling 
bonds (DB) localized on the carbon atoms on the SiC sur- 
face. We notice a splitting between the majority (f) and 
minority spin (!) DB states. In Fig. [2] (b) we show the 
band structure for a FS-12-zGNR, where the convention 
of Ref. is adopted. In Fig. [2{c) we present the dis- 
persion relation for a system with a trench 20. QA wide. 
In this situation, the free-standing carbon atoms form 
a region equivalent to a FS-12-zGNR. The edge mag- 



netism present in the FS-zGNR [30| is also present in 
the system with the trench. One can find both an anti- 
ferromagnetic (AF) configuration, and a ferromagnetic 
(FM) one. Similarly to the free standing nanoribbons. 
we also find that the AF has a lower energy than the FM 
configuration, and the energy differences decay quickly 
with the width of the ribbon, again as obtained in FS- 
zGNR. For the AF case, we compare the energy bands ol 
a FS-12-zGNR with the dispersion relation for the sys- 
tem with the trench (shown in the Fig. [SJb-c)). The 
main characteristics, i.e., the presence of a gap of simi- 
lar value, and the shape of the bands around the Fermi 
level, are very close in both systems. The DB levels that 
are present in the system without the trench still appear 
even in the presence of a trench, since their origin comes 
from the non- free-standing region. In Fig. HJd-c), the 
corresponding bands for the same systems, but when a 
ferromagnetic configuration is adopted, are shown. One 
can note that as it occurs for the FS-zGNR, the system 
with a trench becomes metallic. 

The geometry and the local magnetization for the AF 
case are presented in Fig. [TJb-c) - side and the top 
views, respectively - where one can identify DB states 
with a non-zero local magnetization as well as local mag- 
netic moments on the the free-standing graphene atoms 
located just above the wall of the trench, which mimic 
the edge atoms when compared to FS-zGNRs. An analo- 
gous behavior can be noted in the FM configuration (Fig. 
[TJd-e)), however with a reversal of the magnetization on 
one of the edges. Thus, similarly to what is observed 
in FS-zGNRs, the magnetic ordering (AF or FM) is de- 
fined by the local magnetic moments that appear at the 
graphene edges, above the wall of the trench, precisely 
the sites that mimic the edge sites of FS-zGNRs. More- 
over, the local magnetization coming from the DB states, 
that appear between the substrate and the graphene, 
can be flipped up or down without changing the semi- 
conductor or the metallic character of the AF and FM 
states, respectively. However, we find a small magnetic 
coupling between the free-standing region edge sites and 
their nearest DB, that favors an anti- ferromagnetic con- 
figuration (if one DB is flipped an increase of O.OieV in 
the total energy occurs). 

It is well know that in systems with localized states, 
which is the case investigated in this work[31], the self- 



interaction correction [331 ] is important to obtain a bet- 
ter quantitative description of the energy bands. Thus, 
to include such correction, we implement the method- 
ology proposed by Filippetti et al. 34] in the SIESTA 
code[35|. In Fig. [3] (f) we show the energy bands in- 
cluding the self-interaction correction for the graphene 
on a SiC substrate without a trench. Comparing with 
the same system without this correction, (shown in Fig. 
[2ja)), one can see that the splitting between the t and 
I DB states increases from 1.5 eV to 2.6 eV. Without 
the self-interaction correction, the highest occupied state 




FIG. 3. (Color Online) Energy bands in the Y — M graphene 
direction, (a) Single sheet of graphene over a crystalline SiC 
surface, (b) Armchair FS-GNR with 7 hexagons from one 
edge to another, (c) Graphene layer over SiC containing a 
trench 18.2 A wide along the Y — M graphene direction, (d- 
f) Self-Interaction-Correction calculaions for the same cases 
presented in (a-c). The solid (blue) lines correspond to the t 
states, and the dashed (red) lines are for the 4- states . 



is a DB state, whereas after the inclusion of such cor- 
rection it becomes a substrate bulk state. Other effects 
of this correction include the increase of the local mag- 
netization at the edge atoms of FS-zGNRs (from 0.28 
to 0.37 Hb), and a considerable increase in the gap. In 
Fig. [2fh) we present the energy bands for the system 
with a trench in the T — K direction with the inclusion of 
the self-interaction correction for the AF configuration. 
This calculation can be well compared with an AF FS- 
12-zGNR, shown in Fig. d](g), where the self-interaction 
correction is also employed. 

If we now consider the trench along the T—M graphene 
direction, the system mimics an armchair FS-GNR, as 
can be seen in the Fig. [3] For the sake of compari- 
son we plot the dispersion relation for a single layer of 
graphene over a SiC surface without the presence of a 
trench along the T — M direction. In Fig. EJb) we show 
the band structure for a FS-15-aGNR (following the no- 
tation of Ref. 24) . In [3fc) we show the energy bands 
when the graphene is over a SiC surface containing a 
trench of width 18.2 A. The flat bands in Figs. [3](a) and 
(c) at energies —0.5 eV and 1.0 eV are due to the DB on 
the surface of the SiC. The energy bands presented in 
Figs. [21b) and (c) are very similar to each other around 
the Fermi energy. The gap for the FS-aGNR is slightly 
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FIG. 4. Energy gaps. In (a) the zig-zag case, and in (b) the 
armchair case. In both graphs calculations using the Self- 
Interaction-Correction are considered. 



smaller, and the graphene layer on SiC presents a small 
spin splitting of the valence band. The energy bands 
shown in Fig. G^d-f) include the self-interaction correc- 
tion, and uses the same geometries of Figs. Efa), (b) and 
(c), respectively. Interestingly, it has no effect on the 
gap; the DB states are essentially moved away from the 
Fermi level in such a way as to restore the spin degen- 
erate system and to bring the behavior of the trenched 
graphene nanoribbon closer to the free-standing case. 

We have also calculated the gap dependence with the 
width of the trenches. In |Ua) we present the results 
along the V — K graphene direction. The free-standing 
regions of these systems are similar to the zigzag FS-iV- 
GNRs, with N = {4, 6, 12, 16}. In this case, the energy 
gap monotonically decreases with the width of the trench, 
in a similar fashion to what occurs in the free-standing 
systems 24J. In Fig. 0|b) we consider the T—M graphene 
direction. Similarly to what happens in the FS-aGNRs, 
3 families of systems, named 3p, 3p+l and 3p + 2, can be 
identified. Since the graphene layer is strongly attached 
to the substrate, once the trench direction is defined in 
the SiC surface, there is no possibility for any further 
rotation of the graphene layer. Thus, the orientation of 
the trench relativelly to graphene uniquely defines the 
properties of the ribbons, and intermediate orientations 
of the trenches will lead to edges that have a mixture of 
zigzag and armchair ribbons. 

Summarizing, a single graphene layer on a trenched 
[0001] SiC surface mimics the energy bands around the 
Fermi level as well as the magnetic properties of FS- 
GNRs. With the trench along the T — K graphene direc- 
tion, we find a semiconductor (metallic) state with an AF 
(FM) order, exactly as it occurs in the FS-zGNRs. With 
the trench along the T — M graphene direction, we find 3 
families of systems, similarly to FS-aGNRs. This behav- 
ior occurs because a single graphene layer over a SiC sur- 
face loses the graphene-like properties, and the presence 
of the trench restores the graphene-like properties only in 
a confined strip region. In this way, we propose a natural 
way to create quantum-confined regions in graphene in a 
system that allows large-scale device integration. 



t 
t 

[1] 

[2 
[3. 
[4 
[5 
[6] 
[7] 
[8 
[9 

[10 

[11 
[12 

[13 
[14 
[15 

[16 
[17] 



iiif.usp.br 



mplimaC 

aj rsilva® if . usp . br 
fazzio@if . usp .br 

A. K. Geim and K. S. Novoselov, Nature Mater. 6, 183 
(2007). 

A. H. Castro Neto et al, Rev. Mod. Phys. 81, 109 (2009). 
Y.-M. Lin et al, Science 327 5966 (2010). 
I. Meric et al. Nature Nanotech. 3 654 (2008). 
X. L. Li et al, Science 319 1229 (2008) 
M. Y. Ham et al. Phys. Rev. Lett. 98, 206805 (2007). 
S. Wang et al, Phys. Rev. Lett. 100 206803 (2008). 
Y.-W. Son et al. , Nature 444, 347-349 (2006). 

A. K. Singh and B. I. Yakobson, Nano Lett. 9, 15401543 
(2009). 

D. Gunlycke et al, Nano Lett. 7, 3608-3611 (2007). 

B. Huang et al., Phys. Chem. C 112 13442-1344 (2008). 

C. Berger et al, J. Phys. Chem. B 108 19912 (2004). 

F. Xia, et al, Nano Lett. 10, 715 (2010). 
Y. Q. Wu et al, Appl. Phys. Lett. 92 092102 (2008). 

G. Gu et al, Appl. Phys. Lett. 90, 253507 (2007). 
P. Hohenberg, W. Kohn ,Phys. Rev. 136, B864 (1964); 
W. Kohn and L.J. Sham, Phys. Rev. 140, A133 (1965). 
We use the SIESTA codefTl] with the Troullier-Martins 
pseudopotential parametrization 19j], a mesh cuttoff 
of 170 Ry for the real space grid integration, and 
a k-sampling grid of 1 x 1 x 7, where the 7 k- 
points are taken along the trench direction. The GGA 
exchange-correlation functional is employed with the 
PBE parametrization (2(J. All geommetries were fully re- 
laxed within a force criterion of 0.03eVyA, and a vacum 
interval of 15 A between the slabs is used to avoid unde- 
sirable interactions between the periodic images. 

[18] E. Artacho et al, Phys. Status Solidi B 215, 809 (1999). 
[19] N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993- 
2006 (1991). 

[20] J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. 

Lett. 77, 3865 (1996). 
[21] Y. Chen, T. Kimoto, Y. Takeuchi, R. K. Malhan and H. 

Matsunami, Jpn. J. Appl. Phys. 44 4909 (2005). 
[22] P. W. Waldrab, J. W. Kretchmer, J. D. Galea, U. S. 

Patent No. 0117188 (2010). 
[23] S. Shivaraman et al., Nano Lett. 9, 3100 (2009). 
[24] Y.-M. Son, M. L. Cohen, and S. G. Louie, Phys. Rev. 

Lett, 97, 216803 (2006). 
[25] S. Sonde et al, Phys. Rev. B 80, 241406(R) (2009). 
[26] A. Mattausch and O. Pankratov, Phys. Rev. Lett. 99, 

076802 (2007). 

[27] F. Varchon et al, Phys. Rev. Lett. 99, 126805 (2007). 

[28] J. Hass et a/.,Phys. Rev. B 75 214109 (2008). 

[29] J. Hass et al., Phys. Rev. B 78 205424 (2008). 

[30] T. B. Martins et al, Phys. Rev. Lett. 98, 196803 (2007). 

[31] The self-interaction correction changes the Hubbard 

Parameter U of the graphene, and with effect, 

changes the band structure and the value of the edge 

magnetization [3^ ] . 
[32] M. P. Lima, A. J. R. da Silva and A. Fazzio, Phys. Rev. 

B 81, 045430 (2010). 



5 



[33] J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 
(1981). 

[34] A. Filippetti and N. Spaldin, Phys. Rev. B 67 125109 
(2003). 

[35] The implementation of the self-interaction correction in 
the SIESTA code was done by M. P. Lima, A. J. R. da 



Silva, and A. Fazzio, in the same spirit of Pemmaraju et 

al. [36j. This implementation has been successfully tested 

and applied in other works [13]. 
[36] C. D. Pemmaraju et al, Phys. Rev. B 75 045101 (2007). 
[37] J. T. Arantes et al, J. Phys. Chem. B, 113, 5376-5380 

(2009). 



